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1.0  INIftODOCriON 


Host  atmospheric  mljili  routines  mad  numerical  models  require 
watar  npor  information  to  ba  inpat  at  speeifie  ataoapharic  1 ovals. 
However,  whan  watar  vapor  data  ara  ainiaal  it  nay  ba  necessary  to  rely 
on  total  praeipitabla  watar  valnaa  as  ths  only  noi stare  information 
available.  Snob  a  condition  is  especially  likely  to  ocenr  over  the 
oceans,  where  all  upper-air  data  ara  soaree.  The  analyst  than  faces  the 
problem  of  assigning  valnaa  of  mixing  ratio  (or  soma  related  quantity) 
to  specific  levels  of  the  atmosphere  baaed  on  an  approximation  to  the 
integral  of  water  vapor  over  all  levels. 

Several  papers  have  diaenased  the  relationship  between  preeipitable 
water  and  surface  daw  point  temperature.  Smith  (1966)  suggested  that 
the  mixing  ratio  profile  can  be  represented  by  the  power  law 
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where  w  is  mixing  ratio,  p  is  pressure.  X  is  an  empirically  derived 
constant,  and  the  subscript  o  refers  to  the  surface.  Bquation  1  ean  be 

used  to  relate  preeipitable  water  to  w  at  all  levels.  While  X  is  - 

l  For 

computed  fox  a  given  latitude  and  season,  the  solution  for  w  is  Jtl 

constrained  to  a  power  law  distribution.  ..*ea 
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Spooled. 


2.0  EIGENVECTOR  APPROACH 

As  alternative  aetkod  aay  be  weed  to  relate  total  veter  to  airing 
ratio  at  apeeifie  level a.  Tbia  aetbod  la  baaed  on  eoapntiag  tbe 
eapirical  orthogonal  fanotiona  (aiganvaotora)  of  airing  ratio  pcofilaa. 
Oivan  a  data  baaa  of  profiler  for  a  particular  latitnda  band  and  aaaaon, 
aiganvaotora  can  be  eoapntad  which  repreaent  the  baaie  aodaa  of  vertioal 
variation  in  water  vapor.  The  firat  eigenvector  rapraaanta  the  pr Leary 
node,  which  daacribaa  how  varlationa  in  total  water  tranalate  into 
airing  ratio  profile  variationa.  Kendall  (1973)  givaa  a  ooaplete 
diacnaaion  of  aiganvaotora  and  the  aaaoeiatad  prinoipal  ooaponanta. 
Vachtaann  (1973)  givaa  an  araapla  of  the  naa  of  aiganvaotora  to 
repreaent  ataoapheric  profiler. 

Any  airing  ratio  profile  veotor  w  nan  be  rapr a ranted  aa  a  linear 
coabination  of  the  aiganvaotora  and  a  veotor  £  of  ooaff icienta. 

1  “  E«  +  &  <*> 


where  Jg  ia  a  square  aatrir  with  aaeh  oolaan  coapoeed  of  one  eigenvector, 
and  w  is  the  aaan  airing  ratio  profile.  If  we  ignore  all  ainor  aodaa 
of  airing  ratio  variation,  we  can  approriaate  £  wring  only  the  firat 
eigenvector  a  and  ita  correaponding  ooaffioiant  a 
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Proa  Bqwation  3  it  ia  olaar  that  aa  approriaation  to  the  airing 


ratio  profile  aaa  be  obtained  by  finding  the  appropriate  eoeffieieat  a. 


give  a  that  the  firet  eigeaveotor  tad  the  aeta  profile  are  epeeified  ia 
adwaaee . 

The  reaaiaiag  task  ie  to  ehow  how  e  eaa  be  deteraiaed  f roa 


preeipitable  water  0.  Preoipitable  water  ie  defiaed  ae 


If  the  trapesoidal  approxiaatioa  ia  weed  for  the  iategral  thea  this 
beeoaee 


1L- 

« - ;  *  v*i 

•i-i 


where  Wj  -  (wj  ♦  w^jJ/lj  ip^  “  Fj  **  F|_j»  ui  L  ie  the 
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at  whieh  aeaeareaeata  are  aade.  Sabatitatiag  froa  Bqaatioa  3  giwee 
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eiace  a  -  Z  v  .Ap.  aad  b  ■  Z  e.Ap.  u«  eoaetaate  for  a  givaa  latitude 
i-1  ai  1  i-1  1  1 

band  aad  aaaaoa.  The  vat ar  vapor  profile  oaa  tbaa  ba  eaaily 
approximated  froa  preoipitable  vatar  bp  aaiag  Equation  t  aad  then 
Equation  3. 

3.0  RESULTS 

The  mixing  ratio  aaaaa  aad  the  firat  eigeaveetora  vara  ooapatad  at 
the  aix  aaadatory  lavala  for  five  doaaiaa  of  aaaaoa  aad  latitude  aaiag  a 
data  bare  of  ooaaaio  radioaoada  aeaaureaenta  of  aixiag  ratio.  Table  1 
defiaaa  the  five  doaaiaa  aad  givaa  the  ooapatad  valuer  of  a  aad  b.  The 
firat  eigeaveetora  are  givaa  ia  Table  2  aad  the  aeaa  aixiag  ratioa  are 
givea  ia  Table  3. * •* 

Table  1.  The  five  doaaiaa  for  ooapatatioa. 


Naaber  ia 


Doaaia** 


Latitade 


.-1  .-2 


-  TM 

30°S-30°lN 

All 

173 

680.2 

4 

-  MLS 

30°-60°N 

Saaaier 

184 

508.1 

4 

-  MLV 

30°-60°k 

Viator 

35 

325.9 

3 

-  POS 

«0°-90°k 

Saaaer 

97 

255.7 

4 

-  POW 

60°-90°K 

Viator 

116 

137 .8 

3 

*Suaaer  ■  1  Nay  -  31  October,  winter  ■  1  Noveaber  -  30  April 

•*TR  ■  Tropioal,  ML  -  Mid-latitude,  PO  -  Polar 


Doaaia 


Table  2.  Firat  Eigeaveetora. 
Preaaure  Level  (kPa) 


Explaiaed 
Variaaee  (%) 
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Table  3.  Mean  Ml zinc  Eatios  (g/kg) 


Pressure  Level  (kPa) 


Domain 

30 

40 

50 

70 

IS 

100 

1-1* 

.22 

.63 

1.45 

3.96 

9.0S 

13 .91 

2  -  MLS 

.19 

.59 

1.17 

3.18 

6.30 

3  -  MLV 

.15 

.47 

.87 

1.77 

6.97 

4  -  POS 

.09 

.36 

.76 

1.91 

3.04 

4.39 

5  -  POI 

.02 

.13 

.33 

.97 

1.75 

2.56 

The  aean  nixing  ratios  and  the  eigenvectors  can  be  coapared  using 
the  graphs  in  Figures  1  and  2,  respectively.  As  expected,  there  is  acre 
total  water  in  summer  than  in  winter,  and  aore  in  lower  latitudes.  For 
each  doaain  the  first  eigenvector  shows  the  vertical  distribution  of 
mixing  ratio  variance  within  the  first  aode  of  water  vapor  variability. 

The  influence a  of  the  aeans  and  eigenvectors  for  the  five  doaains 
are  deaoaatrated  by  the  ourvea  in  Figure  3.  The  prof ilea  shown  were 
derived  by  solving  Equations  t  and  3  for  an  arbitrary  value  of  0  »  2.22 
g/ea2,  and  uaing  the  appropriate  values  of  e,  j^,  a,  b  for  each  doaain. 
Note  that  the  ahapea  of  the  polar  auaaer  and  winter  (POS  and  POP) 
prof ilea  are  uarealiatie  in  thia  case  beeauae  2.22  a  is  a  very  high 
value  of  preeipitable  water  for  the  polar  region.  The  profile  shapes 
are  actually  a  different  function  of  preeipitable  water  for  eaoh  doaain. 
4.0  COMPOTE*  1E0GKAM 

A  FOMftAN  computer  program  has  been  written  which  performs  the 
conversion  of  preeipitable  water  values  into  approximate  aixing  ratio 
profilea.  All  of  the  neceaaary  eoeffieienta  are  stored  in  the  program, 
so  the  only  needed  inputs  are  latitude,  Julian  date,  and  preeipitable 
water.  The  code  ia  brief  and  includes  ample  documentation.  The 
Appendix  to  this  report  provides  a  program  liating  and  sample  output. 


Figmre  2.  Foe  auk  of  the  fir*  Itulai  listed  is  Table  1,  the 
plot  shove  the  first  eigesreotor  of  the  aisiag  ratio 
oovsrisaee  aatrix. 
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PROGRAM  PHP*OFL( IN PUT (OUTPUT (TAPES- INPUT# TAPE6-QUTPUT l 
******  THIS  PROGRAM  DERIVES  A 

******  »hoST  PROBABLE"  MIXIN6  RATIO  PROFILE  FRQN  A  PREC1PITA6LE  HATER 
******  VALUE.  THE  COEFFICIENTS  USEO  ARE  SPECIFIC  TO 
******  LATITUDE  BANO  AND  TINF  OF  YEAR. 

VARIABLE  IDENTIFICATION 

******  e  •  COEFICIFNTS  OF  THE  FIRST  EIGENVECTOR 
******  TOP  ■  NUMERATOR  COEFFICIENT 
••••••  OEN  ■  DENOMENATOR  COEFFICIENT 

******  VM  •  MEAN  MIXING  PATIO  PROFILES 
******  v  ■  DERIVED  NIXING  PATIO  PROFILE 
*»♦»**  PH  •  PRECIPITABLE  HATER 

******  ISESN  ■  SEASON  INDICATOR  (HIMTER-1*  SUMMER— 1) 

♦*♦*»*  I  •  IDENTIFIES  SEASON/LATITUDE-BAND  DOMAIN 
C  •  FIRST  PRINCIPLE  COMPONENT 


**  STORAGE  ALLOCATION 
DIMENSION  TnPI5)*0EN(5)#HN(6*3)*V(6l*E(6#9l*P(6l 
DATA  TOP/680. 2*508. 1,329. 9*255. 7*137. 8/ (DEN/42. 90*41. 90* 36. 08# 

♦  49. 67# 38. 34/# P/30. *40. *90. *70. *85. >100./ 

OATA  UM/OO. 22* 00. 63*01. 45*03. 96*09. 08* 13. 91* 

♦  00.19*00.99*01.17*03.18*06.30*10.47* 

♦  00. 19*00. 47*00. 07*01.77*04. 03*06. 97* 

♦  00.09*00.36*00.76*01.91*03.04*04.39* 

♦  00.02*00.13*00.33*00.97*01.75*02.56/ 

DATA  E/.011* .038*. 099* .302# .945# .773# 

♦  .009* .031*. 083* .285, .973**763* 

♦  .009* .036*. 073*. 127*. 940*. 828* 

♦  .003# .059# .148# .433* .670* .982# 

♦  .000*. 013. .093*. 244*. 929* .813/ 


*****  WRITE  A  HEADER  ♦» 

WRITE (6# 080 1 (P(LI#  L-l#6) 

080  FORMAT! "1*#»LATITU0E  JULIAN  PHICH)”#10X*”HIXING  RATIO  (6/KGI* 
♦  /13X»"DAY»,13X*6(F6.0#2X)#"«— PRESSURE  (KPAI-) 


*♦***•  PEAO  A  OATA  SET  AND  IOENTIFY  CLASS 
DO  420  J-l*500 
READ! 9 *110) PLAT. JDAY* PH 
110  FORM«T(1X,F6.2*I10*F10.2) 

IF (EOF ( 5) .NF .OIGOTO  440 
ISESN— 1 

IF  UJOAY.LT.  121). OP.  (JDAY.GT. 304  mSESN«l 
IF (PLAT .GE.O. IGOTO  210 
PLAT..PLAT 
ISESN— ISESN 

210  IF(PLAT.6T. 30. IGOTO  220 

I-l 

GOTO  290 

220  IF(PLAT.GT. 60. IGOTO  240 
IF( ISESN. GE.O)  GOTO  230 
1-2 

GOTO  290 
230  1-3 

GOTO  290 

240  IF( ISESN. GE .0)  GOTO  250 
1-4 

GOTO  290 
250  1-5 
290  CONTINUE 


*****  COMPUTE  MIXING  RATIOS  AND  PRINT 
C-(196.2*PH-T0*(II I /DENI  I) 

00  320  L-l*6 

W(L)-HM(L*TI«C«E(L*Ii 
320  CONTINUE 

HPITE(6#410)  PLAT* J DAY# PH# (Hill#  L-l#6) 
410  F0RMAT(1X*F6.2*U0,F10.2*6F8.2) 

420  CONTINUE 


440  STOP 
ENO 


hi  ^  m  ^  m  9  j*  9  m  ^  m  m  -  • 


.*  %  *•  %  "i  *p  »*•  «  •  .  '  j*  ‘  i  ■  •  •  •  *•*'■'* 


ff.f. 


